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Abstract 

We calculate in the SU(6) gauge theory the mass of the lightest flux loop that winds 
around a spatial torus, as a function of the torus size, taking care to achieve control of 
the main systematic errors. For comparison we perform a similar calculation in SU(4). We 
demonstrate approximate linear confinement and show that the leading correction is consistent 
with what one expects if the flux tube behaves like a simple bosonic string at long distances. 
We obtain similar but less accurate results for stable (fc-)strings in higher representations. 
We find some evidence that for k > 1 the length scale at which the bosonic string correction 
becomes dominant increases as N increases. We perform all these calculations not just for 
long strings, up to about 2.5'fm' in length, but also for shorter strings, down to the minimum 
length, l c = 1/T C , where T c is the deconfining temperature. We find that the mass of the 
ground-state string, at all length scales, is not very far from the simple Nambu-Goto string 
theory prediction, and that the fit improves as N increases from iV = 4 to N = 6. We estimate 
the mass of the first excited string and find that it also follows the Nambu-Goto prediction, 
albeit more qualitatively. We comment upon the significance of these results for the string 
description of SU(iV) gauge theories in the limit N = oo. 



1 Introduction 



There is considerable numerical evidence for linear confinement in SU(2) and SU(3) gauge 
theories and significant evidence that the effective string theory governing the long-distance 
dynamics of confining flux tubes is in the same universality class as the simplest Nambu-Goto 
bosonic string. (See [UI2IS] for some recent calculations.) For SU(iV > 3) gauge theories, while 
there is some evidence for linear confinement [HHj, these questions have not been addressed 
with much precision. In particular, this is so for the stable fc-strings in higher representations 
which have been the object of intensive recent study (HIHEIIH!- One usually assumes a simple 
bosonic string correction for these /c-strings, and one also assumes that this leading correction 
dominates on the same length scales as it does for strings in the fundamental representation 
in SU(2) and SU(3). Since the comparison with theoretical predictions for fc-string tensions 
is sensitive to this assumption, it is important to check its validity. In addition not much is 
known about the effective string theory at large N. Since it is in the N = oo limit that the 
gauge theory might be equivalent to a string theory [2j and since it is in that limit that one 
can currently make quantitative contact between field theories and their string theory duals 
jH], it is important to learn what are the properties of the effective confining string theory in 
that limit. 

In this paper we address these questions with a lattice calculation in the SU(6) gauge 
theory. From other work on the mass spectrum and the deconfining temperature [HJ EH IE]] 
we know that N = 6 is very close to N = oo for many quantities. This is no surprise since 
the leading large-iV correction is expected to be 0{1/N 2 ). Our calculations are at a fixed 
value of the lattice spacing a. In units of the fundamental string tension a, the value of a is 
a ~ 1/4-y/er which is well into the weak coupling region where any corrections to continuum 
physics should be small. For comparison we also perform a calculation in the SU(4) gauge 
theory at a similar value of a. 

The main questions we address are as follows. 

• Do we have linear confinement at large iV? To address this we calculate the energy of a 
'flux tube' that winds around a spatial torus. We then vary the length of that torus and see 
if the energy grows (approximately) linearly with the length, once the length is large. 

• What is the leading correction to this linear behaviour at large /? Equivalently: what is the 
central charge of the effective string theory that describes the long distance confining physics? 

• If we reduce the length of the flux tube - indeed, to its minimum value - how does its energy 
compare to the behaviour predicted by the simple Nambu-Goto action? This addresses the 
question of what is the effective string theory on all scales. In particular, do things become 
simpler as N — > oo? 

• How does the energy of the first excitation of the string compare to the Nambu-Goto 
prediction? This probes the content of the effective string theory in a different way. 

• How does all this affect the reliability of our calculations of fc-string tensions? And indeed, 
what is our best estimate of these tensions? 

In the next Section we shall describe in greater detail the background to these questions. 
We then summarise the technical details of the lattice calculation. We then move on to the 
calculation itself. We begin by quantifying what are potentially the most serious systematic 
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errors in such calculations. We then interpret our results. We finish with a summary of what 
we have learned and then describe the next step in improving upon these calculations. 

2 Preliminaries 

We shall begin with a discussion of how our calculations of the string energy as a function of 
length can probe the nature of the effective string theory. We then briefly outline the lattice 
setup of the calculations. 

2.1 Strings 

If we have linear confinement, we expect that the mass of a flux loop that winds around a 
spatial torus of length I will be given by [TTl E] 

m(l) l =-al-^ (1) 

where c is the central charge of the effective string theory describing the long-distance prop- 
erties of the confining flux tube. (The dimensional factor, d — 2, has been included in the 
coefficient.) There is significant numerical evidence that c = 1 for SU(2) and SU(3) gauge 
theories PUIS El ; corresponding to the universality class of the simplest bosonic string. In this 
paper we shall attempt to determine the value of c for SU(6), in the expectation that this will 
also be its value at N = oo. 

There are old ideas that in its confining phase the SU(iV = oo) gauge theory might be a 
string theory. It is therefore interesting to calculate mil) not just for large I but for all /, and 
indeed for excited strings as well. This is also interesting in the context of the more recent 
dual string approaches: while c in eqn(£Q) determines the number of exactly massless modes in 
the string theory, there will be massive modes generated in the descent to a boundary SU(oo) 
gauge theory, which may be encoded, for example, in the excited state spectrum of the string. 

The very simplest bosonic string theory in d space-time dimensions, based on the Nambu- 
Goto action, has an energy spectrum [J3J E] 

£ „ W = ff( { 1 + ^(„_^)}' ( 2 ) 

and a ground state energy, in d=4, 

The correction to the universal 0(1/1) term in E (l) is 0(1/P) and, as one would expect, 
the expansion is in inverse powers of l 2 a. This is consistent with the observation in PQ that 
the universal 0(1//) term is the dominant correction for ly/o > 3. If one adds the most 
general higher- dimensional operators to the free term in the action, while remaining within 
an effective string theory with the same degrees of freedom, then one can show [T3| that the 
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expansion in eqn(j2J) retains its form except that the coefficient of the 0(l/l 3 ) term is now 
undetermined. If one includes in addition to interactions amongst the massless transverse 
modes, massive degrees of freedom which naturally arise in dual string constructions, then 
these would typically be encoded in excited states whose energies do not approach E (l) 
as I — > oo jinj. Thus we can expect the string spectrum to contain interesting and useful 
information about the details of the equivalent string theory. In this spirit we will calculate 
the energy of the lightest and first excited strings as a function of 

Of course, as one decreases / one encounters a minimum length for which a winding flux 
loop exists. This is easy to see. Suppose we choose the flux loop to wind around the x-torus. 
Now imagine relabelling co-ordinates x,y,z,t — ► t,x,y,z. In this case the length I equals 
the inverse temperature and the mass of the winding string gives the effective string tension 
(times I) [Jj]]. If T c is the deconfining temperature then for T > T c there is no confinement 
and no winding string. Thus the minimum value of I is l c — 1/T C . 

A linearly confining theory naturally embodies deconfinement because the number of 
strings of length I grows as oc exp{+7/} with 7 > and this will clearly overwhelm the 
Boltzman suppression oc exp{— crl/T} once T is large enough jTTj. This argument naturally 
suggests a vanishing effective string tension at T = T c s i.e. a second order transition. As we 
see from eqn(J2J), this occurs in the Nambu-Goto action at the length / = 1/T C S where E (l) 
vanishes, i.e. at 

y/o ~ y (d-2)7T ~ [0.977 : d = 3 U 

Of course this string deconfining temperature might be trumped by a transition to the gluon 
plasma at a lower value of T = T c < T c s , driven by the large entropy of that plasma. In this 
case it is natural for the transition to be first order. It is therefore interesting to compare the 
predicted values in eqn(HJ) to the calculated values of T c in those cases where the transition is 
in fact second order. One finds [TBI EE] 



T 

- 1 - C 



0.709(4) 

1.12(1) 

0.98(2) 



SU{2),d = A 

SU(2),d = 3 (5) 
SU(3),d = 3 



We note a quite remarkable co-incidence with the Nambu-Goto values. (This remark is not 
new: see e.g. jTHj.) Of course what we are really interested in is the N — ► 00 limit, where 
there are theoretical arguments that the gauge theory might be equivalent to a string theory, 
and here deconfinement is known to be first order jTH]. Nonetheless this provides us with a 
bound 

lim -f- > lim — - = 0.596(4) d = 4. (6) 

The Nambu-Goto string is of course extremely simple - no doubt much too simple. From 
string and dual string approaches [T3J one expects 

T s 1 

-7- « -r 7 
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where c counts the number of effective degrees of freedom in the string theory, on the scale 
of T c . Thus one typically would expect a value of T^/^/a that is lower than the Nambu-Goto 
value. And if one could obtain limiv-»oo T^/^/a one might learn something interesting about 
the dual string theory. It might seem that this limit is inaccessible given that the transition is 
first order for all N > 3. However this is not necessarily so [20] • If the tension of the interface 
between confined and deconfined phases is oc N 2 (its natural dependence) then the hysteresis 
around T = T c can become large E] on even moderately large spatial volumes. That 
is to say, if we are in the confined phase we will stay in that phase as we increase T above 
T c and even as T — > T c s , if the hysteresis is extended enough. This is a realistic possibility. 
Preliminary calculations [20] in SU(8) indicate that for N = 8 the hysteresis is not (yet?) 
wide enough. Nonetheless the calculation provides an improved bound of T^/^/a > 0.63. 

One can extend the above discussion from strings in the fundamental representation to 
strings in higher representations pQ. It is useful to classify strings by how their sources 
transform under a gauge transformation belonging to the centre, z G for SU(iV), since 
gluon screening does not transform sources between these classes. Let the transformation be 
z k . Then the lightest string in this class is called the fc-string. (The usual string is thus a 1- 
string.) One can think of it as a collection of k fundamental strings and an interesting question 
is whether they form a bound string state whose tension, cr^, is less than the k separate strings, 
i.e. <jfc < ka. It is also easy to see that k is less than the integer part of N/2 so that we have 
to go to at least SU(4) to have the possibility of this new type of string. In fact this possibility 
is realised pQ: the /c-strings are strongly bound and the spectrum has been calculated quite 
accurately for k < 4; N < 8 [6|. Since the Nambu-Goto string action has no interactions, it 
cannot give such a binding and so cannot be the whole story at finite N. However as iV —>■ oo 
the binding disappears, — > ka, and the Nambu-Goto action re-emerges as a possibility. 
One might also entertain the possibility of the deconfining transition occurring in discrete 
steps at different temperatures, with the fc-string deconfining at T = T c fc , but it turns out that 
in practice this does not occur - there is a single first order transition for iV > 4 [21]. And 
one can also argue that the same is the case for the would-be string deconfining temperatures 



The calculation of requires calculating the mass of a A;-string of length / and then 
applying a correction as in eqn(Q). Typically one assumes that this correction is the same as 
for the fundamental string, c = 1, and that it is the only important correction for l^Jo > 3 
because that is what one finds in SU(2) and SU(3) for the fundamental string. This assumption 
is important at the level of accuracy required to distinguish between competing theoretical 
ideas about the fc-strings. Thus one of the things we will do is to try to obtain some direct 
numerical evidence with accurate calculations in SU(6) for the /-dependence of the masses of 
the lightest k = 1, 2, 3 strings. 

2.2 Lattice 

We work on periodic hypercubic L x x L y x L z x L t lattices with lattice spacing a. The degrees 
of freedom are SU(iV) matrices, Ui, defined on the links I of the lattice. The partition function 
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Z(J3f) = ! ]JdU l e-^ 1 '^ Tl ^ ; f5 
J i 

where u p is the ordered product of matrices around the boundary of the elementary square 
(plaquette) labelled by p and g 2 is the bare coupling. This is the standard Wilson plaquette 
action and since the theory is asymptotically free and since the bare coupling is a running 
coupling on length scale a, the continuum limit is approached by tuning (3 = 2N/g 2 (a) — > oo. 
One expects that for large N the value of a is fixed in physical units (e.g. in units of the 
mass gap) if one keeps the 't Hooft coupling A(a) = g 2 (a)N fixed i.e. (3 oc iV 2 . This has been 
confirmed |Hj in non-perturbative lattice calculations. 

We will consider a loop of flux that winds around the x-torus, so that it is of length / = aL x . 
(Our notation is that / denotes a length in either physical or lattice units, depending on the 
context, while the capitalized form is always in lattice units.) A generic operator (pi that 
couples to such a periodic flux loop is an ordered product of link matrices along a space-like 
curve that winds once around the x-torus. Correlations are taken in the t direction so that 
the energy of the loop is an eigenstate of the Hamiltonian (transfer matrix) defined on the 
xyz space. Such a correlator may be expanded 

C(t = ant) = <<#M0)> = E l^l0(O)|^ac)| 2 e- a ^ ; E, < E i+1 (9) 

71=0 

where, if we have confinement, E is the energy of the lightest flux loop. Since the fluctu- 
ations that determine the error in the Monte Carlo calculation of C(t) may be expressed as 
a correlation function with a disconnected piece, the error is approximately independent of t 
while C(t) itself decreases exponentially with t. Thus one needs operators 4>i that have very 
large projections onto the desired state(s) so that this state dominates C{t) at small t. This 
can be achieved using now standard techniques (see e.g. [B]). If we have linear confinement E 
becomes large for long loops, and especially so for higher representations, and C(t) becomes 
very small. It can then be difficult, even at small t, to make the error on C(t) much smaller 
than its value, if one uses a standard Monte Carlo. In some of these cases we therefore use a 
recently developed error reduction technique [22] • 

We choose L y = L z for convenience. We need to make sure that L y is large enough for 
there to be no significant finite volume effects. We will also need to choose L t large enough 
so that the only state with a significant probability to propagate around the t-torus is the 
vacuum state. In practice it turns out that for large L x one can use L y = L z = L t = L x , but 
for small L x we shall need to use L y , L t 3> L x . Controlling these different finite volume effects 
is one way in which the present calculation improves upon previous calculations. 

Our calculations will be performed for a single lattice spacing. In the case of SU(6) we 
choose (3 = 25.05. This corresponds to 

a(D = 25.05) = °J^!A = <hl^B. : SU( 6) (10) 
in units of the fundamental string tension a (see below) or the deconfining temperature T c 
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16 . In SU(4) we choose (3 = 10.90, which corresponds to 



The value of the lattice spacing has been chosen to be (almost) the same for SU(4) and SU(6) 
so that we can separate the iV-dependence from the a-dependence (to leading order). This 
value of a is small enough that we can expect lattice corrections to be negligible for our 
purposes. 



3 Results 

We calculate the string mass as a function of its length on a variety of lattice volumes. In 
Table H] we list the masses obtained in SU(6) at /3 = 25.05 for the lightest k = 1,2,3 strings 
and, on the largest lattices, for the first excited k = 1 string. In Table El we list a less extensive 
set of results obtained in SU(4) at (3 — 10.9. 

We will begin by describing the important sources of systematic error in such calculations 
and we show how the results listed in Table ^ enable us to control these errors. Having 
obtained the string energies Eq(1) for I ranging from close to its minimum possible value 
upwards, we address the questions posed in the Introduction. 

3.1 Systematic errors 

There are a large number of potential systematic errors. Here we focus on three that we 
believe are potentially the most important ones in this kind of calculation. 

3.1.1 excited state contamination 

Our variational calculation produces operators $f (t) that have very good but not perfect over- 
laps onto the ground string states. We therefore calculate the correlator C(t) = (<3>^(t)$f (0)) 
and extract the energy from the larger values of t = ant where the excited state contamina- 
tion has died out. The signal that we have reached such t is that the correlator is given by a 
single exponential. A problem arises if Eft (I) is large - either because I is large or because k 
is large. In this case the correlator will be very small beyond the smallest values of t, so that 
it becomes dominated by the statistical error and one cannot obtain any significant evidence 
that the higher excited states have in fact died away at the t at which one extracts the energy. 
Clearly the bias will be to extract energies at lower t as the energy increases and so to risk 
a larger admixture of even heavier excited states, and so to overestimate the energy as k or 
I increases. This risk is enhanced for larger k because the best overlap decreases roughly as 
the /c'th power of the k — 1 overlap. (Which provides a motivation for our description of the 
fc-string as a loosely bound state of k fundamental strings.) It is also enhanced for large I by 
two facts. First the gap between ground and excited states decreases as oc l/l and second 
the overlap onto the lightest string appears to decrease exponentially with its length I. The 
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second effect is not unexpected and is in practice not severe: the overlap reduces from ~ 0.985 
for I = 8 to ~ 0.960 for I = 16 using a very similar basis of operators. 

All this is not an issue for the smallest values of / where we have accurate calculations of 
C(t) over a large range of t. To obtain some control of this error for larger Z, in particular 
for I = 16 and Z = 20, we have performed not only a normal Monte Carlo calculation of the 
correlator, but also a separate calculation using a novel and powerful error reduction technique 
|22[ I23j. The latter calculations are starred in Tabled The error reduction technique works 
best for a particular mass range for which it needs to be tuned in advance. Here we have 
tuned the calculation to work best for masses close to those of the lightest k = 2 string. The 
unstarred calculations use a standard Monte Carlo algorithm and, although the masses are 
extracted from an exponential fit to C(t) between t = a and t = 4a, in the case of Z = 16 and 
Z = 20 the signal-to-error ratio grows so rapidly with increasing t that, effectively, the masses 
are obtained from t — a to t — 2a. The starred calculations involve extracting the masses 
from t = 2a to t = 3a where we can be quite confident that any contamination from heavier 
excited string states will be negligible. 

Comparing the results from the starred and unstarred calculations in Tabled we find very 
good agreement within the pairs of calculations for all the states. This provides good direct 
evidence that any contamination from excited states is negligible in the present calculation. 

3.1.2 transverse size corrections 

Calculations of winding flux loops have usually been performed on L 3 L t lattices. Is a transverse 
size of L large enough for a flux tube of length L or does this create significant finite volume 
corrections? To the extent that L is large compared to the width of the flux tube, and the 
oscillations are simple single-valued harmonic modes, one would expect the corrections to be 
very small. So unlike the systematic error discussed above, this is expected to affect short 
rather than long strings. And to the extent that a fc-string is like k weakly bound k = 1 
strings, we might expect any error to be more severe for larger k. 

To address this question we start by comparing the three calculations in Table Q that 
have L x = 8 and L t = 30. (The corrections induced by too small a value L t will be discussed 
below.) These three calculations have transverse sizes L y = L x = 8, L y = 10 and L y = 12. We 
observe that there is no discrepancy in the masses obtained with the last two, implying that 
for L x = 8 a transverse size of L y = 10 is in fact large enough. The L y = L x = 8 calculation, 
on the other hand, produces masses that are very different: while for k = 1 the shift is a 
modest ~ 10%, for k — 2,3 the shifts are ~ 24% and ~ 35% respectively. Our expectation 
that the corrections should be larger for larger k thus appears to be confirmed. 

Moving to L = 10 we have L y = 10 and L y = 12 to compare. For k = 1 there appears to 
be no difference, while there is evidence for a modest discrepancy ~ 5% for the k = 2 string. 
That this is not merely a statistical fluctuation is confirmed if we compare the k = 2 effective 
masses obtained from the correlator between t = and t = a, where the relative statistical 
error is much smaller. The conclusion is that for L = 10 there is a shift in using L y = 10 but 
that it is very small. Moving on to L = 12 we compare L y = 12 and L y = 14. The lightest 
and first excited k = 1 strings show no difference and the possible difference for the k — 2, 3 



7 



strings is shown to be no more than a statistical fluctuation when we look at the statistically 
much more accurate t — 0, a effective masses. 

We therefore conclude that, for this value of a, once L > 12 it is safe to use an L 3 L t lattice 
for such string calculations. Using a^/a ~ 0.25 this translates into the statement that if the 
length I = aL x satisfies l^/o > 3 then an equal transverse size is adequate. It so happens 
that, for quite different reasons, lattices used for calculating £;-string tensions have satisfied 
precisely this criterion. 

3.1.3 temporal size corrections 

If all the spatial dimensions, Lj, are large then all the string energies will be large, and 
exp{— EL t } will be negligibly small for all states other than the vacuum if we choose, for 
example, Lf — Li. (Recall that all glueball masses will in practice be large compared to 1/L t .) 
So for larger Lj the partition function is indeed dominated by the vacuum and there is no 
problem. For small L, however, we can only be sure of excluding the spatially periodic strings 
from propagating right around the temporal torus if we choose L t so that exp{— E (L)L t } <C 1. 
This means that we must increase L t as we decrease L. If we do not do so, then the partition 
function, which is a normalising factor within all our calculated expectation values, will no 
longer dominated by the vacuum contibution and we run the risk that our various calculated 
energies will no longer be with respect to the correct vacuum energy. 

In practice one finds that for the k — 1 string such an effect is quite weak - as we see by 
comparing the 8 3 14 and 8 3 30 lattices in Table^or the 10 4 with the 10 3 16. This is as expected: 
e.g. for the 8 3 14 lattice one sees that exp{— Eo(L)L t } ~ 0.015 and there is a factor of 3 for the 
strings in the different spatial directions. So the expected shift in the effective vacuum energy 
is only AE vac ~ 0.003. In fact it is well known that even with smaller L t the corrections are 
smaller than expected, and this can be ascribed to a cancellation between these extra modes 
in the path integral numerator and in normalising denominator (partition function) in the 
calculation of an expectation value. One would expect such a cancellation if the interaction 
between this extra mode and the string propagating between the operators of the correlator 
was weak. 

However, as we see in Tabled there is a quite dramatic mass shift for the k = 2 and k = 3 
masses when we compare, for example, the 8 3 14 and 8 3 30 lattices. We can interpret this, in 
part, as follows. Consider our k = 2 calculation. We have two k = 2 operators separated by t 
and we assume that for large enough t the correlator is dominated by the lightest k = 2 loop 
propagating between these operators. However an alternative contribution arises from two 
k = 1 loops propagating between the operators in opposite directions around the temporal 
torus. Such a contribution is t-independent and acts, in the calculation, like an unexpected 
non-zero vacuum expectation value for the k = 2 loop operator. This will lower the calculated 
effective mass at larger t. It is therefore clear from this that we require L t to be large enough 
that exp{— E^ =1 (L)Lt} <C exp{— Eq (L)t} for any value of t at which we might evaluate the 
mass of the fc-string. 
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3.2 Interpretation 

Our above discussion of systematic errors suggests that the following choice of lattices from 
TableQwill provide reliable string masses. For I = 7, 8, 10 we use the 7 x 16 2 x 40, 8 x 12 2 x 30 
and , 10 x 12 2 x 16 lattices respectively. For / > 12 we average over the pairs of calculations 
at each I. 

In FigQ we plot the lightest and first excited string masses versus their length /. The 
first thing we observe is that the lightest loop mass, E (l), grows more-or-less linearly with 
its length, /. But before claiming that this provides direct evidence for linear confinement we 
need to establish that this is not some sub-asymptotic behaviour at short distances. Now, if 
we extract a string tension (see below) we obtain a^/a ~ 0.25. Thus our largest flux loop 
length is 20a ~ ~ 2.5'fm'. (The 'fermi' units are obtained using the QCD value of the 

string tension and are only introduced to provide an intuitively familiar scale of length.) This 
is long enough to make it plausible that we are indeed seeing the asymptotic linear rise of a 
linearly confining SU(6) gauge theory. Since we expect corrections to the SU(oo) theory to 
be 0(l/iV 2 ) and since the extracted value of the string tension is very similar to that in, for 
example, SU(3) when calculated in units of, for example, the mass gap E|, this provides 
very plausible evidence for linear confinement in the N = oo gauge theory. 

We now attempt to fit the flux loop mass with a leading linear term and a 0(1//) string 
correction as in eqn([TJ. Since these are supposed to be the leading terms at large I, it is not 
surprising that we do not obtain an acceptable fit over the whole range of I. We therefore 
drop lower I values from the fit and find that we first obtain an acceptable x 2 P er degree of 
freedom, x%i f° r I — 10, with a slightly better fit for I > 12: 

a £o (0=A/-^=/ a06383(38), - 1 16(8W3 ' (12) 
3/ [0.06356(47)/- 1.09(10)tt/3/ : I > 12 , x% = 0.70 

Defining an effective value, c e ff, of the string coefficient obtained between neighbouring values 
of/ 

vrrl 1| aE (l 2 ) aE (h) 

Ceff 3U'qJ-—T 2 /T~' ( } 

we show in Fig. El how it varies with increasing /. All this provides significant evidence that 
the effective string theory for long flux tubes is in the universality class of a simple bosonic 
string theory, c = 1. 

There will of course be higher order corrections in 1// and, under certain natural assump- 
tions, the next correction will be down by l// 2 [Hj. If we fix the coefficient of the 1// term to 
the bosonic string value 7r/3, then we find that we can obtain an acceptable fit to our whole 
range of /: 

aE (l) = 0.06358(20)/ -J- : j > 7 , X % = 0.9 (14) 

This is remarkable since, as pointed out earlier, there is a minimum length for a periodic 
flux loop, which is l m in = l/aT c ~ 6.66 in the present calculation. Thus the fit in eqn (fT4^) 
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essentially works all the way down to the shortest possible strings. (If we fit the 0(1/1) term 
as well, then its coefficient comes to 0.94(16), with a slightly worse value of X%-) Note that the 
natural dimensionless expansion parameter is 1/al 2 and so we would expect the l/l 3 term in 
eqn (|14jl to have a coefficient that is larger by a factor 0(1/ a) ~ 16 than that of the leading l/l 
correction. This is indeed consistent with what we find. Note also that the 0(l/l 3 ) coefficient 
is about twice as large as the Nambu-Goto coefficient in eqn©. 

In Fig ^ we also plot the prediction of the Nambu-Goto string theory for E (l), as in 
eqnQ with n = 0. This fit has only one free parameter, a, and visually follows the calculated 
values quite closely although, with \% — 3-0, it is clearly not an acceptable fit. However, for 
the comparable SU(4) calculation (discussed below) we find x% — 7-5. This improvement in 
goodness of fit with increasing N does leave open the intriguing possibility that at iV = oo 
the effective string theory might be governed by the simple Nambu-Goto action. 

We also plot the prediction from eqn((2J) for the first excited state, which now has no free 
parameter at all. We see that this tracks the calculated values quite well, except at the very 
smallest I. It is interesting to note that the gap between the energy of the lightest and first 
excited string state is large, and for most values of / larger than the mass of the lightest 
glueball, which has a mass am,Q ~ 0.74. For such I the excited string can presumably decay 
into the lightest string and a glueball. This should show up in the correlator through an 
effective energy that gradually decreases with increasing t. The fact that we see no sign of 
this is presumably due to the large- N suppression of decays. For SU(4), where decays are not 
so suppressed, we do in fact find it difficult to extract excited string energies for intermediate 
I, and it is possible that this is part of the reason. 

In FigJS] we plot the masses of the lightest k = 2 and k = 3 strings. These are heavier 
than the k = 1 string and thus possess larger errors. Nonetheless the approximate linear 
growth with / is evident. We plot fits with the bosonic string correction and we list the best 
coefficients of a l/l correction in Table El The k = 3 calculation is not accurate enough 
to be informative on this issue, but for the k = 2 string we have some statistically weak 
evidence that for I > 12a ~ 3/a/o" the dominant correction is 0(1/1), with a coefficient that 
is consistent with that of the simple bosonic string. To investigate this in more detail we 
analyse the k = 2 string masses using eqn ()13)) . and plot the result in Fig. [21 Comparing the 
values of c e // from the k = 1 and k = 2 strings, we see that the k = 2 values are much larger 
for shorter strings, and, while consistent with approaching the bosonic string value, probably 
does so for somewhat larger values of I. This behaviour is, in fact, not unexpected if we think 
of the k = 2 string as two loosely bound k = 1 strings. If it acts as two incoherent, unbound 
strings it receives twice the string correction that it would receive if it behaved as a coherent, 
single string. Since the correction is negative and grows with decreasing I, it is clear that 
there should be some length /' such that for / < /' it will be energetically favourable for the 
fc-string to be k incoherent k = 1 strings rather than a single coherent A;-string. Moreover 
since — > ko as N — > oo, it is also clear that I' — > oo as iV — > oo jSJ. More generally we can 
think of the Liischer correction as inducing a repulsive interaction between the two strings, 
which leads to a value of c e // that smoothly decreases from 2 to 1 at a length scale that grows 
with N. It is tempting to see Fig. 121 as displaying such a behaviour. 

In addition to the above SU(6) calculation we have also performed, for comparison, a less 
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extensive and less accurate SU(4) calculation, as listed in Table El A fit with a 0(1/1) bosonic 
string correction plus an extra 0(1/ 1 3 ) correction gives 

aE (l) = 0.05798(32)/ - ^ - 11 2 ^ : /> 8 ,$ = 1.0 (15) 

As with SU(6) the coefficient of the 0(l/l 3 ) is of the expected order of magnitude, although it 
is no longer possible to include the smallest value of I in the fit. In Fig0]we plot the lightest 
and first excited k = 1 flux loop masses and compare to the Nambu-Goto string predictions. 
As we remarked earlier, the best Nambu-Goto fit is considerably worse than for SU(6). It is 
thus possible that such a fit will become acceptable in the N = oo limit. 

Fitting the k = 2 string masses with an 0(1/1) bosonic string correction and an additional 
0(1 /I 3 ) term, we obtain as acceptable fits, 

aE k=2 (l) = f0.07963(53)/ - vr/3/ + 1.7(21.2)/Z 3 : I > 10 , SU(4) 
" l! [0.10739(75)/ -tt/3/-65.7(10.6)// 3 : l>8 , SU(Q) 

The fact that the 0(1 /I 3 ) coefficient is much larger for iV = 6 than for N = 4 presumably 
encodes the fact that the minimum string length at which the Luscher correction becomes a 
good approximation is larger for larger N. (And the fact that it is much larger than for k = 1 
in eqn((TU) suggests it grows with k.) Together with eqns (|14H5|) . these lead to 

Vkaeo.1t/T, f 1.373(12) : SU(A) 

o " \ 1.689(13) : SU(Q) [ ' 

These ratios have very small systematic and statistical errors and show quite clearly that at 
this fixed but small value of the lattice spacing these string tensions lie in between the Casimir 
scaling [21] and 'MQCD' conjectures [25] with which they are usually compared [UHOE]. 

As an aside, we note that the above calculations also tell us about how the effective string 
tension, <J e ff, varies with temperature T. Let us relabel our Euclidean axes so that our x 
direction is the t direction. Then our k — 1 Wilson line operator is just the world line of a 
static fundamental source and if we define a e ff to be the coefficient of the linear piece of the 
free energy of two such static sources then it is clear that [TBJ [21] 

l=l/T Ep(l) l=l/T 7T 2 4 

a eff(T) = —j— = °-t^ t +0(T) (18) 

if the string is in the simplest bosonic universality class. In Fig]S]we plot the values of a e ff(T) 
against T/T c for both SU(4) and SU(6). There is a hint that as iV increases the effective 
string tension has smaller corrections at all T. 

4 Conclusions 

In our calculation in the SU(6) gauge theory we have obtained good evidence for linear con- 
finement of flux loops that are up to / ~ 5/y/a long. This is long enough in physical units (it 
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corresponds to ~ 2.5fm in QCD) for us to see it as good evidence for exact linear confinement. 
The calculation is at a fixed lattice spacing but this is small enough, a ~ that we are 

confident that we are seeing continuum physics. Finally iV = 6 is large enough that we see 
this is as good evidence for linear confinement in the N — > oo limit. 

For strings in the fundamental representation our SU(6) calculation shows that the leading 
0(1/1) correction at large I has the coefficient one expects if the effective string theory belongs 
to the simplest bosonic string universality class. This coefficient is determined accurately 
enough for this result to be significant. 

These calculations were performed with care to control and quantify the main system- 
atic errors. We found that string tensions could be reliably calculated on L 4 lattices once 
L > This fortunately coincides with the conventional choice of lattice size for such 

calculations lj. For flux loops in higher representations our results do not pin down the 
string corrections very accurately, but they are consistent with a bosonic string correction 
dominating at large I. There is, however, evidence that this occurs at larger I than for the 
k = 1 string, the more so the larger N, and that the effective string coefficient at smaller 
I is larger for larger k. This fits in with the picture of a k string as being k loosely bound 
fundamental strings with the Liischer correction inducing an effective repulsive interaction. 
It also implies that past calculations of £;-string tensions may have suffered from a small but 
systematic downward bias. The ratio of the k = 2 and k = 1 string tensions that we obtain in 
this paper lies more-or-less mid-way between the Casimir Scaling and the 'MQCD' motivated 
trigonometric formula. 

Remarkably, we found that using only the 0(1/1) bosonic string correction with an ad- 
ditional 0(1 /I 3 ) term, we could describe the (fundamental) string energy all the way down 
to the minimum string length, I ~ 1/T C where T c is the deconfining temperature. The fitted 
coefficient of the 0(l/l 3 ) term turns out to be 0(l/a 2 a) just as one would expect on general 
grounds. 

The string energy is also found to be close to the prediction of the simplest string theory, 
that governed by the Nambu-Goto action, and the comparison between SU(4) and SU(6) 
leaves open the possibility that for SU(oo) the effective string theory is indeed governed by 
that action - although this would be hard to understand on theoretical grounds EJ. Our 
calculation of the energy of the first excited string and our discussion of string deconfinement, 
reinforces the conclusion that whatever effective string action governs the N = oo theory, it 
must be 'quite close' to Nambu-Goto. 

It would be useful to have, say, SU(4) and SU(8) calculations that are as precise as our 
SU(6) ones, so that we could directly extrapolate the string energy to iV = oo as a function of I, 
all the way down down to the minimum length, I ~ 1/T C , dictated by 'spatial deconfinement'. 
In fact, as we remarked earlier, one can use the hysteresis of this first order phase transition to 
go to lower I. The aim would be to go all the way down to the second order string condensation 
transition where the string mass vanishes. Although it appears that the hysteresis is not wide 
enough for this to be possible in SU(8), there is good reason to believe it will be possible for 
not very much larger values of N. The low-/ behaviour of the string energy together with the 
excited string spectrum, particularly for states whose energies do not approach the ground 
state energy at large I, should provide useful insights into the detailed nature of the effective 
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string theory for, and the dual string theory to, the SU(oo) gauge theory 
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Figure 1: The masses of the lightest, •, and first excited, o, k — 1 flux loops that wind around 
a spatial torus of length I in the SU(6) calculation at /3 — 25.05. The dotted lines are the 
predictions of the Nambu-Goto string action, as in eqn(J2J). The dynamical lower bound on 
the string length is l min = l/aT c ~ 6.63. 
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Figure 2: The effective coefficient of the l/l universal string correction, calculated from 
eqn ()13|) . for the range of lengths indicated. For k — 1, •, and k — 2, o, strings in SU(6). 
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Figure 3: The masses of the lightest k — 2, •, and fc = 3, o, flux loops that wind around a 
spatial torus of length I in the SU(6) calculation at (3 — 25.05. The dotted lines are the best 
fits with a bosonic string correction, as in eqn(JTJ) with c = 1. 
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Figure 4: The masses of the lightest, •, and first excited, o, k — 1 flux loops that wind around 
a spatial torus of length I in the SU(4) calculation at (3 — 10.9. The solid line is the best fit 
with a bosonic string correction, as in eqnQ with c — 1. The dotted lines are the predictions 
of the Nambu-Goto string action, as in eqn(J2J). The dynamical lower bound on the string 
length is l min = l/aT c ~ 6.66. 
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Figure 5: The finite temperature effective string tension a e ff(T) plotted as a function of T/T c 
for SU(6), •, and SU(4), o, at a fixed lattice spacing a ~ 0.15/T C . 
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